Preparation of Nonclassical States in Cavities with a Moving Mirror 
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We describe how a quantum system composed of a cavity field interacting with a movable mirror 
can be utilized to generate a large variety of nonclassical states of both the cavity field and the 
mirror. First we consider state preparation of the cavity field. The system dynamics will prepare a 
single mode of the cavity field in a multicomponent Schrodinger cat state, in a similar manner to 
that in a Kerr medium. In addition, when two or more cavity modes interact with the mirror, they 
may be prepared in an entangled state which may be regarded as a multimode generalisation of even 
and odd coherent states. We show also that near-number states of a single mode may be prepared 
by performing a measurement of the position of the mirror. Secondly we consider state preparation 
of the mirror, and show that this macroscopic object may be placed in a Scrodinger cat state by a 
quadrature measurement of the light field. In addition we examine the effect of the damping of the 
motion of the mirror on the field states inside the cavity, and compare this with the effect of cavity 
field damping. 
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I. INTRODUCTION 

Recently, Mancini et.al. B have shown how a cav- 
ity with a movable mirror (treated as a quantum har- 
monic oscillator) can be used to synthesize Schrodinger 
cat states (||| of the cavity field. In fact, as we will 
see here, this system can lead to the production of an 
extensive class of novel states of the cavity field, includ- 
ing entangled states of two or more cavity modes. One 
of the aims of this paper is to describe how these can 
be achieved. Furthermore, in Ref. W only the effect of 
photon leakage from the cavity was considered as a rel- 
evant source of decoherence. We examine the opposite 
extreme, namely when the decoherence induced in the 
field due to the damping of the mirror's motion is the 
dominant source of decoherence and photon leakage is 
almost absent. We calculate the effect of this type of de- 
coherence on the states of the cavity field and explicitly 
demonstrate that, due to the oscillatory nature of the sys- 
tem, the decoherence rate is much lower than expected. 
Moreover, we show that apart from trying to improve 
the mirror's isolation, increasing its frequency also helps 
to reduce the rate of the type of decoherence considered 
here. This fact comes as a benefit in the generation of 
atleast one type of nonclassical state of the cavity field. 
These results definitely brighten the prospects of observ- 
ing nonclassical states of the field inside the cavity given 
that damping of the mirror's motion is inevitable. It is 
also shown that even when the effect of mirror damping 
on the cavity field is dominant, it does not destroy all the 
nonclassical features of the cavity field. 

In Ref. 0], the main focus was the generation of non- 
classical states of the cavity field. We give here a thor- 
ough treatment of the quantum dynamics of the mirror 



as well. In particular, we point out ways in which this 
mirror motion can also be put in a nonclassical state. Of 
course, if these are to be observed, then very good iso- 
lation of our system from the environment is necessary. 
However, the very fact that this system allows in principle 
the production of nonclassical states of a macroscopic ob- 
ject such as the mirror (after sufficient isolation) should 
be interesting in itself. 



II. DYNAMICS OF THE UNDAMPED SYSTEM 

Numerous authors have previously treated the system 
of a cavity field and a movable mirror, both quantum 
mechanically [QJ^] and classically and both theoret- 
ically jl],(|,[§ and experimentally |6|7||. We here assume 
the movable mirror to be a quantum harmonic oscillator 
with frequency Lu m and annihilation operator denoted by 
b, interacting with a cavity field mode of frequency ujq 
and annihilation operator denoted by a. The relevant 
Hamiltonian H is 

H = hujQ a) a + Tiuj m b^b — fig a)a{b + b^) (1) 

where 
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and L and m are the length of the cavity and mass of the 
movable mirror respectively. 

The time evolution operator corresponding to the 
above Hamiltonian was evaluated by Mancini et al. Oj 
(for completeness we give a proof in Appendix A) and is 
given by 



U(t) = exp \—ira)a i\ exp [ik 2 (a 1 a) 2 (t — sini)] exp ^ka^airjb^ — exp [— ib^bt\ 



(3) 
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where rj — (1 — e~ 4 *), k — g/oj m is the scaled cou- 
pling parameter, r = wo/w m , and t represents a scaled 
time, being the actual time multiplied by u> m . We note 
that values for k of the order of unity are experimen- 
tal/ feasible (ujq ~ 10 16 s -1 ,u/ m ~ 1 Khz, L ~ 1 m, m ~ 
10 mg) ||J|q' We see from the above equation that 
there is now an explicit Kerr- like term in U(t), so that 
physically one might expect the cavity field to have an 
evolution similar to that in a Kerr like nonlinearity. In 
view of the large variety of nonclassical states that can be 
produced by a Kerr medium ]^-]Io[|, this system clearly 
offers prospects for the production of nonclassical states 
of light Indeed, it is known that, like the Kerr 

nonlinearity, this system of a moving-mirror cavity also 
exhibits optical bistability jlH . 

Let us assume that initially both the mirror and the 
cavity field are in coherent states. To see that this is rea- 
sonable we first note that the long-time steady state of 
the cavity mode is the vacuum, and that the most stable 
pointer state for the mirror (being a single harmonic 
oscillator) is a coherent state Jl2| . The steady state there- 
fore consists of the cavity mode in the vacuum state and 
the mirror in a coherent state. The cavity mode can now 
be placed in a non-vacuum coherent state by driving it 
with a coherent input field on a time-scale which is much 
shorter than the time-scale of the mirror motion. Thus 
we write the initial state at time t — as 

|*(0)) = |a)c ® |/3) m (4) 

where |a) c and |/?) m are initial coherent states of the field 
and the mirror respectively. 
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FIG. 1. The linear entropy, S, of the mirror state, in the 
absence of any damping, is plotted here as a function of time, 
and for various values of the scaled coupling parameter k. We 
have taken the initial coherent state amplitude of the cavity 
mode to be a = 2, and that of the mirror to be (5 — 2. The 
entanglement of the mirror with cavity field increases from a 
value of zero to a maximum at t = n and subsequently falls 
back to zero at t — 2-n. Both the scaled time t and the linear 
entropy S are dimensionless quantities. 



The time evolution of the system in the interaction 
picture (that is, omitting the free evolution of the field) 
leads to a state at time t given by 

2 OC n 

|tt(*)> = e-^^^ e * fc2 » 2 ( i - sW )| Jl ) c8 |^( t )) m (5 ) 
„ \>n\ 

n— 

where |n) c denotes a Fock state of the cavity field with 
eigenvalue n, and the \<j> n {t))m. are coherent states of the 
mirror given by 

\Mt))m = I P^ lt + kn(l - e- u ) ) m . (6) 

Two features of the mirror dynamics emerge from the 
above equations: 

(1) After a time t — 2tt the mirror returns to its orig- 
inal state. At all times between t = and t — 2w the 
mirror state is entangled with the field state with the 
entanglement being maximum when t = it. The time 
dependence of the entanglement is shown in Fig.(l) as a 
plot of the purity, or "linear entropy" [T^ | which is given 
by (1 — Tr(pM(t) 2 )) where pmH) is the reduced density 
matrix of the mirror. 

(2) The mirror, always being in a mixture of coherent 
states during its evolution, is described by entirely posi- 
tive Wigner functions and therefore has a fully classical 
dynamics as illustrated in Fig. (2). In this figure, and all 
the figures showing the Wigner function of the mirror, x 
shall stand for b + and y for — i(b — tf). The mirror 
undergoes an oscillation of a different amplitude (but at 
the same frequency) for each number state of the field. 
The net effect is a kind of breathing of the mirror state; 
the mirror state undergoes an oscillatory increase and de- 
crease of its position and momentum spread. However, 
this is very different from the well known breathing of 
squeezed states, as the mirror state is a mixture of co- 
herent states. 
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FIG. 2. The Wigner function, y) , of the mirror at var- 

ious times t of the system in the absence of damping. We have 
taken the initial coherent state amplitude of the cavity mode 
to be a — 2, and that of the mirror to be (3 — 2. The scaled 
coupling parameter has been set at k — 0.5. The quadratures 
x and y and the Wigner function W(x,y) axe given here in 
dimensionless form. 



Thus, during undamped evolution, no nonclassical 
state of the mirror is generated. However, because of the 
entanglement of the mirror with the field, it is possible to 
generate nonclassical states of the mirror by performing 
conditional measurements on the field. We discuss this 
in section IV. 



III. GENERATION OF NONCLASSICAL STATES 
OF THE CAVITY FIELD 



We divide the discussion of the generation of non- 
classical states of the cavity field into three categories: 
(A) Multicomponent cats generated at t = 2ir due to 
dynamics alone (i.e without any external intervention). 
Among these, a method of generating the 2-component 
Schrodinger cat has been discussed in Ref. fjj]. (B) En- 
tangled states of two or more light modes generated at 
t — 2-k due to the dynamics alone. (C) Nonclassical states 
produced due to conditional measurements on the mirror. 
Of course the discussion here is by no means exhaustive 
and there remains the possibility of generating even more 
interesting states. 



A. Multicomponent cats 



At time t — 2ir the state of the cavity field, as obtained 
from Eq.(5), is given by 



Depending on the value of the parameter k, the state |£) (b) k=0.4082 

can be made equivalent to a variety of multicomponent 
cats. For k = 0.5, 

0.06-. 




where c\ = 1 and c 2 = -c 3 = (1 + e l7r / 3 )/(2i sin(7r/3)) 
(not normalized). For k = l/(2y/2), 

\U)c = ^-(\a) c - | - a) c ) + \{\ia) c + | - ia) B ) (10) 



which is a four component cat (not normalized). Thus 
by adjusting the ratio of the coupling and the mirror fre- 
quency and thereby varying k one can, in principle, ob- 
tain all these types of cats at time t = 2ir. The Wigner 
functions of the cats produced by different values of k 
are shown below in Fig. (3). In this figure and all figures 
showing the Wigner function of the cavity field x = a+a^ 
and y = —i(a-a^) are dimensionless quantities. Exactly 
the same type of states have been noted before as aris- 
ing in the Kerr medium . It is indeed worth noting 
that creating such states in a cavity such as ours has 
intrinsic advantages because there is an extensive set of 
tomographic methods ||l4|-po|| that can be implemented 
to reconstruct the Wigner function in a cavity. 
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FIG. 3. The Wigner functions of the state of the cavity field 
at t — 2ir are plotted here for various values of the scaled cou- 
pling parameter k. We have taken the initial coherent state 
amplitude of the cavity mode to be a — 2, and that of the 
mirror to be f3 — 2. The quadratures x and y and the Wigner 
function W(x,y) are given here in dimensionless form. 



B. Entangled states of two or more cavity modes 

Two modes of light, each separately interacting with 
the movable mirror, and with no direct coupling between 
them, can end up being in an entangled state at time 
t = 2t: depending on the value of the parameter k. Let 
there be two different modes of light inside the cavity at 
t = 0. For simplicity we assume them to have the same 
frequency and hence same value of the parameter k, but 
mutually orthogonal polarization directions. We assume 
them to be initially fully disentangled and prepared in 
coherent states. Let the initial state of the composite 
system of the mirror and the light be 

|tf(0)) = |oi) e i ® |a 2 ) c2 ® \p) m (11) 

where |ai) c i and |«2)c2 are states of the first and sec- 
ond modes respectively, while |/3) m is the initial state of 
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the mirror. The state of the two light modes evolves at 
t = 2ir to the state (not normalized) 



\Ce) 



00 n n m i 2 irk (n+m) 
n,m=0 



n'.y/m 



|n) c i ® |m) c 



(12) 



E — 

m=0 



m! 



The state can be rewritten as 

Kb 
where 

|£(m,fe)) c2 



e lfc2 ™ 227r |m) cl ® |£(m,fc)>e2 



(13) 



(aie 

n=0 



ifc 47rm\n 



|n) c2 . (14) 



As is evident from the above equation, |£(m, fc)) C 2 are 
states of the same type as the state \Q C given by Eq.(7) 
with an amplitude a,\e lk 47rm . Hence, for certain values 
of k (such as 0.5 or 1/(V6) or 1/(2-^/2)), number states 
of one mode become correlated with multicomponent cat 
states of the other mode. However, if k is such that 
(2k 2 nm) is an integer for all n and to, then the state 
|£(m, fc)) C 2 becomes independent of m and \C,e) becomes 
a disentangled state. Some simplification of Eq.(14) for 
k=0.5 results in the entangled state 



|Ce) = (l + *)l + «l)d® |+02)c2 

+ (1 - i)| - ai) c i ® I - a 2 )c2- 



(15) 



Note that production of entangled states in a Kerr 
medium, when the two modes of light have a direct in- 
teraction between them have been discussed before [ pl| ; 
however this case is different in the sense that the modes 
become entangled because they both individually inter- 
act with the same movable mirror. 

There is no restriction on the number of modes that can 
be entangled using the above procedure (only each mode 
to be entangled has to be resonant with the cavity) . It is 
interesting to demonstrate how eigenstates of the opera- 
tor (a c ia c2 ...a c ^) p where a c i,a c2 , ...,a c j\j are annihilation 
operators of N different modes of light and p is an integer 
can be created in the cavity. As it is possible for resonant 
modes to have frequencies that are integer multiples of 
each other, let us consider the case when there are N dif- 
ferent modes of light in the cavity, with their frequencies 
^ci)W C 2) ....,o;cN being related by w C j = i]j u> c i where rjj 
are integers. Then from Eq.(|2|) it follows that the value of 
their k parameters will be related by kj = rjj k%. Hence 
if the initial state of the cavity modes and the mirror is 



l*(0)> 



K)cl <8> \a 2 )c2 ® ■■■■ ® K)cN 



(16) 



the final (at time t = 2ir ) composite state of the cavity 
modes will be 



E 



i2irfci (jii rii+r/2 n 2 + ...+r]N n N ) 



ni,n2,...,njy— 



N 

n 

3=0 



(17) 



As the term e^ 1 ^ 1 ni+ri2 " 2+ - ■■+''" nN ~> 2 cannot be split up into a product of the form fi(n 1 )f 2 (n 2 )...f N (n N ) 1 for 
arbitrary values of fci, \Ce n ) can, in general, be an entangled state depending on the value of k\. Now consider the 
case when 27rfci 2 = n/p. We have 



(a c ia c2 ...a cN ) p |CE JV ) 



E 



i2irfe 1 2 (»7 1 111+1)2 n 2 + . ••+>)« ™jv) 2 



n 1 ,n 2 ,...,n N =p 



N rij 
If ' 



\ n 3 ~P)cj ■ 



(18) 



Redefining rij in the above equation as rij + p , recasting the equation in terms of the new rij and using the fact that 
ijj are integers, one gets 



(a cl a c2 ...a c ^y\CE N ) = {a 1 a 2 ...a N Ye l ^ + ^+- +r '^ 2 \C En ) . 



(19) 



Hence |Ce n ) are eigenstates of the operator {a c ia c2 ...a c ^) p , generalizations of the well known even and odd coherent 
states which are eigenstates of squares of annihilation operators 0] . 



C. Nonclassical states produced by conditional 
measurements on the mirror 

Because of the entanglement of the mirror state with 
the state of the cavity field during evolution, any mea- 
surement on the mirror will project the cavity field to 
some state which we can determine. We consider here 
measurements of the mirror's position, x = b + w (i.e, 
x represents actual position of the mirror multiplied by 



(y /, 2muj m /h), m being the mass of the mirror), at time 
t = 7r, as the entanglement is most pronounced at that 
time. Such a measurement will project the state of the 
light to 



\ V (x)) c =Nj2^/* k2n2 (x\<f>nmn) 
n=0 



(20) 
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Thus for different values of the mirror position x, differ- 
ent states of the cavity field are produced. Let us choose 
parameters k = I, a = 2 and /? = 2. Then the peak of 
the gaussian function (x\4> n (ir)) lies at x n = An — 4 and 
its width is 6 = 1. Thus the distance between the peaks 
of (x\<f>n(n)) and (x\<f> n+ i(%)) , that is, x„ - x n +i = 4 , is 
greater than the width S. Hence, if the value obtained for 
x obtained as a result of the measurement is near x n , the 
most dominant constituent of the state \n(x)) c is |n) c . In 
this way, we may expect to produce nonclassical states 
which are very close to number states. For example, if 
the value obtained for x is near (which is near x\) the 
cavity field is projected to the state (unnormalizcd) 

117(0)) « e- 4 |0) c - 2|l) c + 2\/2e- 4 |2) c + 0(1CT 7 ) (21) 

which is shown in Fig. (4a). Similarly when x is near 
X5 « 16, we get a state very close to the number state |5) 
as shown in Fig. (4b). When x is inbetween x n and x n+ i 
we get essentially a superposition of states \n) and \n+ 1) 
as illustrated in Fig. (4c) for n — 2 (the state in the fig- 
ure is approximately — 2|l) c + 2v2|2) c , not normalized). 
Fock state superpositions and their nonclassical proper- 
ties have been discussed in |22] . 

(a) x=x. 
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FIG. 4. The Wigner functions of the state to which the 
cavity field is projected for various results, x, of a measure- 
ment of the position of the mirror. We have taken the initial 
coherent state amplitude of the cavity mode to be a = 2, and 
that of the mirror to be (3 = 2. The scaled coupling param- 
eter is k = 1, and the measurment is made at time t — n. 
The quadratures x and y and the Wigner function W(x,y) 
are given here in dimensionless form. 



It is evident from the above discussion that the genera- 
tion of the above type of states solely rely on the narrow- 
ness of the coherent state width in comparison to the 
spatial separation between the peaks of the gaussians 
(x|^> n (7r)) and (x|0„+i(7r)). As this spatial separation 
can always be increased by increasing the value of the pa- 
rameter fc, the proximity of the states produced to Fock 
states and their superpositions can always be improved. 

We note that the first two types of states, that is, the 
cats and the entangled states, rely on the Kerr type term 
in the time evolution and hence they can be generated 
in a Kerr medium as well. In particular, generation of 
cat states in a Kerr medium was suggested some time 
ago However, the generation of the third type of 

states (i.e the number states and their superpositions) do 
not appear to have an analogue in the Kerr medium be- 
cause they are entirely dependent on the entanglement 
of the states of the mirror and the cavity field during 
evolution. 
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IV. GENERATION OF NONCLASSICAL STATES 
OF THE MIRROR FROM CONDITIONAL 
MEASUREMENTS 

The system also offers the opportunity to create quite 
different nonclassical states of the mirror by conditional 
measurements on the x quadrature of the light field. As 
the spatial separation between the coherent state compo- 
nents of the mirror is maximum at t = ir, measurements 
at this instant are likely to produce the most nonclassical 
states of the mirror. A measurement of the x quadrature 



G 



of the light field at this instant of time projects the mirror 
state to 



£ a e iM {x]n)lM7r)) ^ (22) 

n=0 



As the \4> n {^))m are coherent states with a different am- 
plitude for each value of n, |$(x)) m is a superposition of 
spatially separated coherent states and as such an entirely 
nonclassical state of the mirror. Moreover, as the param- 
eter k is increased, the separation between the coherent 
components |</>„(7r)) m increases. So by varying k, one can 
control the macroscopic distinguishability of the states 
involved in the superposition. In Fig. (5) we have plotted 
the Wigner function for the state of the mirror produced 
when k = 1 and the measurement result is x — 0. This 
state is distinctly nonclassical and looks somewhat like 
a Schrodinger cat state with unequal amounts of each 
coherent component. 




FIG. 5. The Wigner function of the state to which the mir- 
ror is projected when a measurement of the x quadrature of 
the cavity field at time t = n gives a value 0. We have taken 
the initial coherent state amplitude of the cavity mode to be 
a — 0.8, and that of the mirror to be /3 — 2. The scaled cou- 
pling parameter is k = 1. The quadratures x and y and the 
Wigner function W(x, y) are given here in dimensionless form. 



The importance of the above procedure stems from 
the observation that there is already a considerable 
amount of literature, both theoretical ^3| and experi- 
mental flpifl, focussed on methods to prepare atoms in 
nonclassical motional states. So it is only natural to ex- 
pect that placing a more massive object in such states 
should be an issue of serious consideration. The abil- 
ity to place macroscopic objects in non-classical states 
may even posses applications. Hollenhorst, for example, 
has shown that placing gravitational wave detectors in 
squeezed states leads to higher sensitivity p5[ . Allthough 
schemes relying on momentum transer to a massive ob- 
ject from several microscopic objects (such as atoms or 
neutrons) have been suggested ^6|, these are definitely 
difficult to implement experimentally Thus the scheme 



described in this section offers a possible direction from 
which the production of nonclassical states of a macro- 
scopic object may be approached. Obviously, these states 
will decohere very rapidly, and their realization may en- 
able some tests of decoherence models 27 1. The three 
main issues which will need to be addressed to bring this 
scheme to a practical level are: (1) How well isolated the 
movable mirror can be made , (2) how to perform to- 
mography of the states produced and (3) how to perform 
instantaneous measurements of the x quadrature of the 
field inside the cavity. We leave the treatment of these 
issues for the future. 



V. THE EFFECT OF ENVIRONMENT INDUCED 
DECOHERENCE OF THE MIRROR ON THE 
STATES OF THE CAVITY FIELD 

In Ref. jjj, the decoherence of the cavity field was 
treated as entirely due to the leaking of light through 
the cavity mirrors. Here, we investigate the complemen- 
tary case i.e. when the photon leakage from the cavity 
is almost absent, but damping of the mirror's motion is 
significant. It is somewhat artificial to assume no leak- 
age of light from the cavity, but in principle the rate of 
mirror damping can be made a few orders of magnitude 
more than the rate of damping due to the leakage of light 
by choosing mirrors of sufficiently high reflectivity. We 
do not try to address issues of practical feasibility in this 
paper, however. 

As the state of the cavity field is entangled with the 
state of the mirror at times between and 27r, it is ex- 
pected that the decoherence of the mirror will induce a 
decoherence of the cavity field. However, its effect on 
the state of the cavity field is expected to be nontrivial. 
During the standard decoherence due to leaking of light 
from the cavity, the coherent states of light form the rel- 
evant pointer basis |TIJ . However, in the case considered 
here, coherent states of the mirror are in one to one cor- 
respondence with number states of the cavity field. So 
if decoherence of the mirror forces it towards a coherent 
state basis, it will induce a decoherence of the cavity field 
towards the number state basis. This decoherence will 
leave some imprint on the state produced in the cavity 
at time t = 2n. 

The equation which governs the decoherence of the 
mirror depends on the way the mirror is coupled to its 
environment. Here, for simplicity, we assume that the 
mirror is amplitude coupled to the environment at zero 
temperature, which implies the master equation |2q] 



dp(t) 
dt 



= ~[H,p(t)] + l(2bp(t)tf - tfbp(t) - p(0b) 

(23) 



where H is the Hamiltonian of Eq. ([I]) . Since we are using 
a scaled time t, 7 in the above equation is the usual damp- 
ing constant (that is, the reciprocal of the dissipation 
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timescale) divided by the mirror frequency ui m . Eq.(|23|) 
has the feature that it singles out coherent states as the 
pointer basis O). It is expected that natural decoher- 
ence of macroscopic objects should force them towards 
a coherent state basis (to conform with classical reality), 
and Eq.(|23|) accomplishes exactly this. So we expect that 
the solution of Eq.^ will atleast give the basic features 
of the effect of mirror decoherence on the cavity states 
correctly. A similar justification has been given for the 
use of the above type of decoherence in modelling quan- 
tum measurement in Ref. [ p9| . A more appropriate model 
for the decoherence of the mirror is quantum Brownian 
motion [^Tj, but it is also more difficult to solve for our 
system of interest. Work is in progress to solve the quan- 
tum Brownian motion for our system numerically while 



here we present the analytic solution for the evolution 
of our system when decoherence is described according 
to Eq.(23). This solution is expected to have most fea- 
tures (at least the localization of the cavity field towards 
a number state basis) similar to the solution in the case of 
quantum Brownian motion because both drive the mirror 
state to a mixture of coherent states. 



The technique we use to solve Eq.(p3|) is to apply the 
unitary evolution and the nonunitary (decohering) evo- 
lution alternately for short instants of time At, and then 
take the limit as At — ► 0. For simplicity, we assume 
the initial coherent state of the mirror to be the vacuum 
state. The solution for this initial condition is evaluated 
in Appendix B, and the result is 



PM=e-H 2 V ^^e« 2 ^ 2 - ro2 )( t - sint )e- D (™' i V)(m|®|0 ri (7,t)){0rn(7^)l (24) 
„=o,m=o v n\m\ 

where 

Mj,t) = -^(l-e^+^), (25) 

are the amplitudes of the coherent states of the mirror and 

, fc 2 (n-m) 2 7r l-e-T* . (e^/ 2 ) 4 - 1) { e -(^+i/^-l) 



The term e -o{n,mn,t) in Eq.(§4|) is responsible for de- 
coherence. Note that as 7 — > 0, we have D(n, m, 7, t) — ► 
and </> ra ( 7 , t) — ► <p n (t) where <j> n (t) is that given by Eq.(||), 
with (3 set equal to zero. In other words, the solution 
given by Eq. (G4) reduces to the undamped solution given 
by Eq.(g|). The fact that the mirror is in a mixed state 
and does not return to a pure state at t = 2~k as in the 
undamped case is illustrated in Fig. (6) in a plot of the 
linear entropy of the mirror's state versus time. 



0)0.5- 




FIG. 6. The linear entropy, S, of the mirror state in the 
presence of damping. This is plotted here as a function of 
time, and for various values of the scaled coupling parameter 
k. We have taken the initial coherent state amplitude of the 
cavity mode to be a — 2, and that of the mirror to be f3 — 0. 
We have chosen the damping constant to be 7 = 1. The en- 
tropy fails to return to zero because of the entanglement of 
the system with the environment. Both the scaled time t and 
the linear entropy S are dimensionless quantities. 
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The Wigner function for the light field at t — 2ir 
for various values of our scaled damping constant 7 and 
k = 0.5 (that is, when a Schrodinger cat is expected in 
the absence of any decoherence) is given in Fig. (7). Two 
features concerning the type of decoherence considered 
here become clear from these figures: 

(a) 7=0.001 



0.1 . 



0.05> 



x" 




-10 -10 



(d)y=1 

0.04 ^ 




-10 -10 



FIG. 7. The Wigner function W(x,y) of the cavity field 
at time t — 2n for various values of the motional damping 
constant, 7, of the mirror. We have taken the initial coherent 
state amplitude of the cavity mode to be a — 2, and that 
of the mirror to be j3 = 2. The scaled coupling parameter is 
k — 0.5. The quadratures x and y and the Wigner function 
W(x,y) are given here in dimensionless form. 



(b) y =0.01 




(c) y=0.1 



(1) Note that from Fig. 7(a), it is evident that even for 
a value of scaled 7 as high as 0.001 (that is, unsealed 7 of 
1), we have almost no decoherence of the Schrodinger cat 
at all. Thus, unless the mirror is quite heavily damped 
it has almost no effect on the states produced inside the 
cavity. A possible cause for this is that the efficiency 
of decoherence depends on the separation between the 
coherent states of the mirror corresponding to different 
number states of the field. This separation is not con- 
stant, but oscillates between zero and a maximum value 
of 2k(n — to). The decoherence process, does not, there- 
fore, get the chance to act with as much efficiency as 
it would have if the mirror was not a harmonic oscilla- 
tor. Speaking more mathematically, Eq.(|26|) implies that 
D(n, to, 7, t) is proportional to k 2 j and not just 7. As in- 
creasing the frequency decreases k, it also decreases the 
rate of decoherence even if 7 (which relates to the iso- 
lation of the system) is held constant. In terms of the 
absolute time t (obtained by dividing the scaled time t, 
which we had been using, by u) m ), the time scale of 
decoherence depends on the frequency as (from Eqs.(^|), 
(p6|) and definition of k) 




T d (X (LU m ) 3 



(27) 



Thus controlling uj m offers an alternate way to control 
the decoherence timescale of our system. 

It is worthwhile to mention that an analogous situ- 
ation has been pointed out in Ref. 1 30 1 in the context 



of the Jaynes Cummings model, where atomic sponta- 
neous emmission has a much weaker effect on revivals 
of atomic inversion than the cavity field damping. In 
that case, however, the explanation was quite different, 
namely, that spontaneous emmission is independent of 







the cavity field intensity, while cavity damping is inten- 
sity dependent. This logic will also apply to our system 
when the intensity of the cavity field is large. The term 
which causes decoherence, D(n, m, 7), is completely in- 
dependent of the intensity of the cavity field. So as this 
intensity increases, the photon damping begins to domi- 
nate over the mirror's motional damping, as far as influ- 
encing the final state in the cavity is concerned. 

(2) The interference peak of the Schrodinger cat, which 
is a primary feature of this cat, is not significantly low- 
ered by the decoherence process (See Figs.7(b)-7(d)), 
allthough there is a phase diffusion of the state. The 
reasoning for this is the simple fact that the decoher- 
ence is towards the number state basis, which preserves 
the photon number distribution. The photon number 
distribution at t — 2tt is same for both the decohered 
state given by Eq.(HJ) and the undecohered state given 
by Eq.(||). This means that one of the primary signa- 
tures of the Schrodinger cat, namely, an oscillating pho- 
ton number distribution (that is, the probability of odd 
photon numbers being zero for an even coherent state 
and vice versa) is the same for the final state produced 
at t = 2-7T, irrespective of whether there was decoherence 
or not. So the decoherence considered here is not acting 
to eliminate all nonclassical properties of the field that 
would be produced without decoherence. It does destroy 
the phase information, but maintains the nonclassicality 
of the number distribution. This is precisely the reason 
why the interference peak of the Schrodinger cat does not 
appear to be destroyed even when scaled 7 is as high as 
unity, as depicted in Fig. 7(d). 

So far we have only considered how the generation of 
the Schrodinger cat is affected in the presence of the mir- 
ror's motional damping. Let us now briefly pause to con- 
sider how the generation of the other nonclassical states 
of the cavity field may be affected. The entangled s tates 
of two or more cavity modes mentioned in Section. [II B 
are generated by a procedure identical to the generation 
of the multicomponent cats and hence should have simi- 
lar patterns of decoherence. On t he other hand, the Fock 
states mentioned in Section. Ill C will not atall be affected 



by the type of decoherence considered here, as it drives 
the cavity field towards a mixture of Fock states. More- 
over, even after their production they will be very stable 
because, neither the Hamiltonian evolution, nor the de- 
cohering evolution destroys a Fock state of the cavity 
field. Ofcourse one will have to remember that in a real- 
istic case the photon leakage will also be present, which 
does destroy Fock states |gB[. Lastly, comes the super- 
positio ns of t wo successive Fock states, also described in 
Section. IIIC. These will, ofcourse, be seriously affected 
as the decoherence process will be destroying the coher- 
ence between successive Fock states. As the generation 
of the nonclassical states of the mirror depend crucially 
on the coherence between the different \n) c ® \4> n (t))m 
components of the system's state, it will also be affected. 
In fact, at the very instant of generation of a nonclas- 
sical state of the mirror of the type given by Eq.(22), 



the coherence between any two of its components like 
|0„(7r)) m and \<j) m (ir)) m , will be e -D(n,m n ,^ of the un _ 
damped value. 

We have already mentioned that the decoherence 
timescale of the system may be increased by increasing 
the frequency of the mirror. But whether that actually 
helps in the generation of any of the nonclassical states 
described in this paper requires further scrutiny. For ex- 
ample, the multicomponent cats are generated only when 
k 2 t reaches a certain value. But D(n,m,j,t) depends 
precisely on k 2 t and not simply on time. Thus increas- 
ing the frequency will ofcourse increase the decoherence 
timescale, but also increase the timescale for the pro- 
duction of the cat states proportionately, so that they 
are decohered by the same amount whenever they are 
produced. For the nonclassical states of the mirror de- 
scribed by Eq.(p2|), the situation is somewhat better, but 
not without its problems. By decreasing k (through w m ), 
one not only decreases the rate of decoherence but also 
decreases the spatial separation between the coherent 
state components involved in the superposition. Hence 
the components in the superposition may become more 
coherent, but they also become less macroscopically dis- 
tinguishable . However, in the generation of the super- 
position of two successive Fock states, decreasing uj m re- 
ally does help. The reason is clear from Eq. (J20|) . The 
components of the superposition are Fock states of the 
cavity field and thereby do not depend in any way on 
the parameter k. The only contribution of k comes in 
the amplitudes, and two successive Fock states can al- 
ways be made to have significant amplitudes whenever 
the measurement outcome x is about halfway between 
the peaks of the gaussians (x|</!) n (7r)) and (x|(/< n +i(7r)), ir- 
respective of the value of k. Thus there is atleast one 
type of nonclassical state whose generation can be aided 
by increasing the mirror frequency. 



VI. CONCLUSION 



In conclusion we would like to stress that the work 
presented here offers new prospects for the observation of 
nonclassical features of light in the cavity. Even when the 
decoherence due to the mirror's motion is maximum, the 
nonclassicality associated with the photon number dis- 
tribution in the cavity is preserved. However, nonclassi- 
cality associated with the phase distribution is hindered. 
But here also, one can be optimistic from the viewpoint 
that the effect is not dominant unless the mirror is very 
heavily damped. We have not addressed the problem of 
detecting the states produced in the cavity because there 
exists an extensive literature on this topic ]l^-pO]| . 

The primary aim of further work must be to solve the 
system when the mirror is damped according to quan- 
tum Brownian motion models |ll|]3l| so that one can set 
up explicit limits on the parameters such as the tem- 
perature and mass of the mirror required to observe the 
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nonclassical states in the cavity in their most undeco- 
hered condition. Nevertheless, we expect that the fea- 
ture of decoherence of the field towards the number state 
basis should also be present in that full solution, since 
any physically sensible decoherence model must tend to 
localize the mirror state towards the coherent state ba- 
sis. In addition there are prospects to utilize the system 
for tests of quantum Brownian motion, which can be ex- 
plored further. Also, we have considered starting initially 
with only a coherent state inside the cavity. One might 
be able to produce more interesting nonclassical states 
by starting with nonclassical states. In fact, all that can 
be investigated in the Kerr medium can be done with 
this system. Even more can be done with the mirror sys- 



tem perhaps becuase of the opportunities of conditional 
measurements. 
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APPENDIX A: DERIVATION OF THE TIME 
EVOLUTION OPERATOR FOR THE 
UNDAMPED SYSTEM 



The time evolution operator is given by 

jj g— ira a^— itb*b-\-ikta ^ a(b-\-b ') 



(Al) 



where t is the time multiplied by u) m , k = g/uj m and 
r = LUo/uj m . We now consider a unitary transformation 
using the operator 



T 



_ -ka , a(b'< -b) 



(A2) 



Note that this is a displacement operator for the mirror 
in which the displacement amplitude has been replaced 
by the number operator for the cavity mode. Using the 
Baker-Cambell-Hausdorf expansion |3^] , the effect of this 
transformation is readily shown to be 



TbT^ = b + ka)a, 
T^T* = fo f + kcJa, 
Ta) aT^ = a) a. 

Using the fact that 



(A3) 
(A4) 
(A5) 

(A6) 



for any function /, unitary operator U , and arbitrary set 
of operators {A^}, the effect of the transformation on the 
time evolution operator is easily calculated to be 



TU(t)T^ 



— ira^ at —ib^ bt—ik 2 (a^ a) 2 



(A7) 



Multiplying on the left by T\ and on the right by T, 
we obtain the following expression for the time evolution 
operator: 

TJ(j-\ —ira" at ik 2 (a^ a) 2 t ka* a(b^ —b) —ib^ bt — ka) a(b^ —b) 

(A8) 

(Note that to obtain this expression we have swapped 
various exponentials which contain commuting argu- 
ments.) To obtain the final form of U(t), we need to 



swap the last two exponential factors in this expression. 
To acheive this we note first that the BCH expansion 
gives 

e" l&t Va(&t - b)]e lbH = a) a(tf e"" - be lt ), (A9) 



and again using Eq.(A6) we obtain 



-ibH -ka"<a(b^-b) ib f b _ -fca + a(f) t e~ it -6e it ) 



Multiplying both sides on the right by e 16+6 we arive at 
the relation required to swap the exponentials, namely 

e -ib 1 b e -ka 1 a(b'<-b) _ e ~fco t a(& f e~" -be") b (All) 

We may now write the expression for U(t) as 



(A10) 



£j ^ — ira^at^ik 2 (a^a) 2 t^ka^a(b^—b) 

_ -ko} a(tf e- u -be u ) -ibH 
X c o , 



(A12) 



and to obtain the final expression given in Eq.(|3|), we 
need only combine the arguments of the 3 rd and 4 th ex- 
ponentials, which is readily acheived with the BCH rela- 
tion fl3l. 



APPENDIX B: SOLUTION TO THE MASTER 
EQUATION FOR THE SYSTEM WHEN THE 
MIRROR'S MOTION IS DAMPED 



We need to solve Eq.(p3|), which contains two parts. 
The solution of the first part 



-dT = -n [H ' pit)] > 



is known to be 



p(t) = U(t)p(Q)U\t) 



(Bl) 



(B2) 
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where, U(t) is given by Eq.(||). On the other hand, the 
second part, 

^ = l(26p(t)6t - tfbp(t) - p(t)tfb) (B3) 

is known to transform |Aj)(Aj| as |2j| 

PiilXMXA - (Ai| Aj)^~ e ~ 7 ^ lAie -7 '/ 2 ) (Aj-e -7 ^ 2 ! 

(B4) 



where, |Aj) and |Aj) are coherent states. 
Let the initial state of our system be 

p(0) = \a)(a\ c ® |0)(Q| m , 



(B5) 



that is, we assume the mirror to be in a vacuum state 
initially (for simplicity). As the damping localizes the 
mirror state to a coherent state basis, and corresponding 
to each number state of the field, the mirror is driven to 
a separate coherent state, at any time t the general form 
of the density matrix will be 

P(t) = ^2 Pnm(t)\n) (m\ c (g) \<j)n(l,t))(<p m (j,t)\ m . 

71 — 0,771— 

(B6) 

Thus, we can split the evolution of the entire system into 
separate evolutions labelled by m and n. In other words, 
we break our problem into that of calculating the evo- 
lution of each of the entities p nm {t) |0 n (7,t))(</> m (7,t)| m 
separately. In any small time step At d uring which evo- 
lution occurs only according to Eq.([B2|), the amplitudes 
of the coherent states of the mirror change according to 



0n(7, t) -> (1 - jAt)0„( 7 , t) + iknAt, 



(B7) 



while if evolution occurs according to Eq.(B3), then 



</. n ( 7 ,i)^„(7,*)(l-7AV2). (B8) 
Hence we can set up the following differential equation 



dt 



7 

-i^nvy,t) +ikn- -z4>n{l,t), (B9) 



the solution for which (with initial condition </> n (7,0) 
0) is given by 



0n(7><) 



ikn 
i + j/2 



(1 



-(i+ 7 /2)t 



(BIO) 



Next let us determine how p nm (t) evolves. Following 
Eq.(B2) it evolves as 



n (t\ — ». n (f\ p ik 2 (n 2 -m 2 )At(l~cost) 



(Bll) 



while Eq.(B3) transforms 

Pnm (^) 



>(l-e- 



I'll III 

= pnm^e-^' 1 '-^^ 2 ^ 1 ' 2 , (B12) 
where the smallness of At and the fact that 



Jm(0„(7,t)0 m (7,i)*) = (from Eq. flB10|) ) have been 
use d. He nce, evolving p ram (0) according to Eq.(Bll) and 
Eq.( B12 ) alternately and taking the limit At — > we 
obtain 



Pnm(t) = P„m(0)e* 

Xe -i £l<K(7,t)-<M7,t)| 2 <i* 



(B13) 



Combining Eq.(|B10D, Eq (|B13| ) and Eq.(|BJ), and evalu- 
ating the integral in Eq.(B13) one can obtain Eqs.(|24|)- 
( p6| ) of section^, which constitutes the complete solution 
for the density matrix of the system at an arbitrary time. 
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